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Abstract 

For a pair of real or complex scattering potentials : M —>■ C (j = 1, 2) with support Ij 
and transfer matrix Mj, the transfer matrix of vi + V 2 is given by the product M 2 M 1 provided 
that Ii lies to the left of l 2 - We explore the prospects of generalizing this composition rule 
for the cases that Ii and I 2 have a small intersection. In particular, we show that if Ii and I 2 
intersect in a hnite closed interval of length i in which both the potentials are analytic, then 
the lowest order correction to the above composition rule is proportional to . This correction 
is of the order of if vi and V 2 are respectively analytic throughout this interval except at 
X = i and X = 0. We use these results to explore the superposition of a pair of unidirectionally 
invisible potentials with overlapping support. 
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1 Introduction 

Transfer matrices have numerous applications in a variety of scattering problems in physics and 


engineering BE]. This is mainly because of their composition property that allows for the deter¬ 


mination of the scattering properties of a complicated system from the contributions of its simpler 
constituents. This is most simply described in the standard one-dimensional potential scattering 


0 . 


Consider a possibly complex-valued scattering potential n : M —)■ C with an asymptotic decay 
rate such that |1] 



( 1 ) 
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and let k he & positive real (wave)number. Then every solution of the Schrodinger equation 


— + v{x)^jJ{x) = k‘^'ip{x), x G M, 


( 2 ) 


satishes 

ip{x) -)■ as x ^ ±cx), (3) 

where Aj- and B± are possibly fc-dependent complex coefficients [3]. The transfer matrix M of the 
potential n is a /c-dependent 2x2 matrix that fulhlls the relation [1] 


A 

B 


+ 

+ 



(4) 


Similarly to the S'-matrix, M encodes the scattering properties of the potential v. Recalling that 
scattering solutions of ([2]) are given in terms of the left/right reflection and transmission amplitudes, 
and T, according to 



^ikx ikx ^ — C)0, 




^—ikx _|_ 


for X — 00 , 
for X ^ 00, 


(5) 


and using Q, (j4]), and ([5]), we can express the entries Mij of M as [5] 

Mn = T- R^K/T, = R^/T, 

1 (0) 

M21 = -R^/T, M22 = 1 /T. 

These, in particular, imply that detM = 1. We also note that ([T]) is a sufficient condition for the 
(global) existence of the Jost solutions, '0fc+ = ^^*4 V'fc- = of ([2]), [B]. 

Now, suppose that v can be written as the sum of a pair of scattering potentials Uj : M —)■ C 
(j = 1,2) with the same asymptotic decay property as n, such that the support^ of vi lies to the 
left of that of V 2 - Using Ij to label the support of Vj, we express this condition by ‘Ji -< / 2 ’EI Under 
this assumption, we can relate M to the transfer matrix Mj of Vj according to 


M = M 2 M 1 . 


(7) 


This is the celebrated ‘composition property’ of the transfer matrix. It is also called the ‘group 
property’, because transfer matrices belong to the matrix group SL{2, C) and ([7]) involves the group 
multiplication for this group [2]. The primary aim of this article is to seek for a generalization of 
()7|) to the cases where R n /2 is a hnite interval. Without loss of generality, we can identify the 
latter with [0, i], where £ is a real parameter signifying the length of R HR, and demand that for all 
Xi G R\l2 and X2 G R, Xi < X2. We denote this relation by ‘R R.' Figured] provides a schematic 
description of this condition. To summarize, we wish to generalize the composition property ([7]) of 
the transfer matrix, which holds for R -< R, to situations where R ^ R. 

The main tool that we employ in our treatment of this problem is a recent observation made 
in Ref. [7] which identihes the transfer matrix with the solution of an initial-value problem for a 

^The support of a function / : K C is the smallest closed interval outside which / vanishes. 

X R’ means that for every xi G R and X2 & R, we have Xi <X2- 
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Figure 1: (Color online) Graphs of |ni(x)| (solid blue curve) and |u 2 (a:^)| (dashed purple curve) for 
a pair of potentials Vi whose support li satisfy Ji ^ l 2 - In particular, Ji fl /2 = [0, i] for some £ > 0. 


time-dependent Schrodinger equation. Let [a:_,x+] be any closed interval containing the support of 
V, T± = kx±, and M(r, r_) be the solution of 


i5rM(r, r_) = M(r_,r_) = 1 , 


( 8 ) 


where r G [r_,r+], 

J^{t) := w{t)}C{t), /C(t) : = 


1 


—e 


2zr 


-1 


w{t) := 


v^r/k) 


(9) 


and 1 is the 2x2 identity matrix. Then, M = M(r_|_, r_). In other words, we can express M as 


r'r+ 




( 10 ) 


where is the ‘time-ordering’ operation. Notice also that for all ti,T 2 G [r_,r+] with ti < T 2 , the 
matrix M(r 2 , Ti) is the transfer matrix of the truncated potential: 


•= d{kx — Ti)9{t 2 — kx)v{x) = 


v{x) for kx&[Ti,T 2 ], 
0 for kx^[Ti,T 2 ], 


( 11 ) 


where 6 stands for the Heaviside step fnnction; 6{t) = 0 if r < 0, and 6{t) = 1 if r > 0. In other 
words, 

M(r2, Ti) = =!^exp I —i [ ^(r)(irl. (12) 


' n 

\-ir 


A direct conseqnence of dH]) and (E]) is that M(r) satishes 


(13) 


where the snperscript ‘T’ stands for the transpose of the corresponding matrix. More generally, we 
have 


M(r2,ri) =l7exp<i / J^(r)^(ir>. 


0T2 


(14) 
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2 Generalized Composition Rule 

Consider the case where Ji I 2 and Ji fl /2 = [0, ^ (See Fig. [H) Because v = vi + V 2 , we can take 
[r_, r+] = Ji U / 2 . Moreover, let 

v_{x) := 9{—x)vi{x), Vo{x) := 9{x)9{i — x)v{x), v+{x) := 9{x — £)v 2 {x), (15) 

and I±fi and M±_o denote the support and transfer matrix of respectively. Then, /_ C [r_, 0], 
h ^ [0, , and /+ C [0, r+]. This shows that J_ -< Jq -< /+. We also have +Vo + v^ = V 1 + V 2 = v. 
Therefore, we can use the composition property of the transfer matrices, namely Eq. ([7]), to establish 

M = M+MoM„. (16) 

Next, we introduce the analog of M(ri, T 2 ) for the potentials Vj with j = 1, 2; let Mj(ri, T 2 ) be 
given by replacing v with vj on the right-hand side of flT^ . Then (IT^ implies 

M„ = Mi(0,r_), Mo = M(£,0), M+ = M 2 (r+,£). (17) 

Substituting these in flTB]) and using the fact that 

Ml = Mi(£, r_) = Mi(£, 0)Mi(0, r_), M 2 = M 2 (r+, 0) = M 2 (r+, e)M 2 {i, 0), 


we obtain 

M = M2S(e)Mi, (18) 

where e := ki and 

S(e) := M 2 (e, 0)-iM(e, 0)Mi(e, 0)“^ (19) 

Introducing 

:= •= (20) 
and using (IT^ . we can express S(e) in the form 

S(e) = exp ^(r)^dr}]^ exp i JJ^(r)(ir}] exp ^(r)^dr}]^. (21) 


Notice that according to (jD]) and fl20|) . J^{t) = + J^(r). 

The matrix S(e) contains all the scattering features of the potential v stemming from the fact 
that the support of Vi and V 2 overlap in the interval [0, £]; for £ = 0, Eq. flTT]) gives S(0) = 1 . In the 
remainder of this section, we study the behavior of S(e) for ‘small’ e by examining its expansion in 
powers of e; 


S(e) — 1 + ^ 

n=l 

gW .= s(”)(0), 


cW 

n\ 


;= aTS(e). 


( 22 ) 


In order to compute the coefficient matrices Sq"'^ we apply the above scheme for the truncated 
potentials Vjr{x) := 9 {t — kx)vj{x) where r G ^,k£]. This amounts to changing e in flT^ into an 
independent parameter r; 

S(r) := M 2 (r, 0)-'M(r, 0)Mi(r, 0)"^ (23) 
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We can evaluate Sq"'^ by repeatedly differentiating this equation with respect to r and setting r = 0. 

Taking the derivative of both sides of (1231) an using this equation together with (jHl) and flTdll in 
the resulting expression, we hnd that S(r) is the solution of the initial-value problem: 

fS(r) = [^i(r),S(r)|^, (24) 

S(0) = 1. (25) 

Here an overdot denotes differentiation with respect to r, and for any pair of 2 x 2 matrix-valued 
functions A and B, we have introduced 

[A(r),B(r)|^ ;= A(r)B(r) - B(r)A(r), A(r) ;= M 2 (r)-^A(r)M 2 (r), (26) 

M2(r) := M2(r, 0) = 5^exp |-i J <^(r')(ir'| . (27) 

Observe that for r = 0 the r-dependent bracket [[tI^ reduces to the usual commutator; 
[[A(0), B(0)|g = [A(0),B(0)]. Furthermore, the following useful identity holds. 

^r[A(r),B(r)|^ = [A(r), B(r)|^ + [A(r), B(r)|^ - i[A(r), J^2(r)]B(r). (28) 

Equations fl2TD - flTTI) imply that if the potentials Vj are analytic functions in [0,£], then the 
series (l2^ converges. Because vi{x) = 0 ioi x > £ and V 2 {x) = 0 for x < 0, this condition requires 
that vi{x) = vi{£) = 0 and hm 2 ;_j.o+n 2 (x) = ^ 2 ( 0 ) = 0, which are rather restrictive. It 

is not difficult to see that the convergence of (l2^ can also be established for the cases where Vi 
and V 2 are analytic in (0,£) but are respectively allowed to be discontinuous a.t x = £ and x = 0. 
In the following we focus our attention to potentials Vj with these properties, and, without loss of 
generality, identify the value of Vi (respectively V 2 ) at x = £ (respectively x = 0) with its left limit as 
X ^ £ (respectively right limit as x —)■ 0.) The same applies for the hrst and higher order derivatives 
of vi (respectively V 2 ) at x = £ (respectively x = 0); we identify them with the corresponding left 
(respectively right) derivatives. 

Next, we list some of the properties of the matrix /C(r) that we will use in the calculation of 


They follow directly from the dehnition of /C(r), i.e., ([9]). 

/Co :=/C(0) = *(72 + CTs, /Co :=/C(0) =-2*(Ji, /Co ;=/C(0) = 4((J3 -/Co), (29) 

/C(r) = e-*"^3/Co /C(r) = e'^'^^/Co /C(r) = e-'^^^/Co (30) 

/C(r)2 = /C^ = 0 , /C(Tf = tl = -41, [( 73 , /C(r)] = */C(r), (31) 

[/C(r),/C(r)] = -4*/C(r), [/C(r),/C(r)] = -4*/C(r), [(73,/C(r)] = */C(r), (32) 

where 0 and 1 respectively denote the 2x2 null and identity matrices, and (7j, with i = 1,2,3, 
stand for the Pauli matricesjfl 

Now, we are in a position to compute So"'\ We begin by noting that in view of and ([25]), 

:= S(0) = 0 . (33) 

^According to (l 3 Tl) and (l 32 l) . for all r e M, the matrices Ki := —A 1 (t)/ 4 , K2 := Ai(r)/8, and K3 := (73/2 are 
generators of the algebra sm(1, 1 ); they satisfy [Ki,K2] = iKs, [K2,K3] = iKi, and [K3,Ki] = —*K2. 
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Next, we use fl20|) . ([26]), and the first equation in ([ST]), to establish 


( 34 ) 


If we differentiate both sides of flM|) and use Eqs. fl28|) and fIMl) to simplify the result, we arrive at 

iS(r) = [^i(r),S(r)|^ + [^i(r), S(r)| (35) 

Setting r = 0 in this equation and employing fl25l) and fl33l) yield 

;= S(0) = 0. (36) 

We can similarly carry out the calculation of S*^^^(r) and S^^^(r). Taking the derivative of fl35|) 
and using (1^ and fl34)) give 

= [^i(r),S(r)|^ + 2[jri(r),S(r)|^+ [jri(r), S(r)|^ - i[jri(r),^(r)]S(r). (37) 

In light of ((20]), fl25l) . ((331), ((Ml), and the hrst equation in ((32|), Eq. fl371) implies 

S® := S(^)(0) = -4i'«;i(0)M;2(0) /Cq. (38) 

To determine S*^"‘^(r), we make use of (1^ to differentiate (1371) . This gives 

= [^/^^(r),S(r)|^ + 3[^i(r),S(r)|^ + 3[^i(r),S(r)|^+ [^i(r), 

-i|[^(r),^(r)]S(r) + 2[^i(r),^(r)]S(r) + ^(^[^i(r),^(r)]S(r)j |, (39) 

where we have also benehtted from (l34l) . In view of (ITHl) . (l25l) . (l32l) . (l33l) . (l34j) . (l36l) . and (l38l) . 
Eq. ((39|) implies 

:= S^'^^(O) = -16m;i(0)m; 2(0) (Ji - 4i [m;i(0)m;2(0) + 3u;i(0)m;2(0)]/Co. (40) 

Substituting fl5^ . fl5B]) . fl55]) . and (140]) in flOOD . we hnd the following more explicit expression for 
S(e). 

S(e) = 1 - y M;i(0)M;2(0)/Coe^ 

-^|4u;i(0) m;2(0) (Ji + i[M;i(0) ^ 2 ( 0 ) + 3tf;i(0) m;2(0)]/Co| + O(e^) 

2i 

= 1-—Wi{€)w2{0)K:oe^ 

-^|4M;i(e) W2{0)(Ji + i [wi{e) ^ 2 ( 0 ) - M;i(e) ^ 2 ( 0 )] /Coj + O(e^) 

= 1 - ^^^i(^)^2(0)/Co£^ 

-^{4A:ni(£) ^ 2 ( 0 ) di + z [ni(£) n'(0) - v[{e) ^ 2 ( 0 )] JCo] t (41) 

where 0{e'^) stands for terms of order d and higher in powers of e, and we have used the Taylor 
series expansion of tci(r) about e for r = 0, namely 


m;i( 0) = M;i(e) + 


00 


n=l 


-1 

w^\e)e^ = wi{e) 


n\ 


r7i(e)e + O(e^), 


(42) 


6 





to establish the second equality. 

Equation fl4ip implies that whenever both Vi{€) and ^ 2 ( 0 ) are nonzero, the leading order correc¬ 
tion to the standard composition rule for transfer matrices is of the order of This correspond to 
situations where vi and V 2 are respectively discontinuous a.t x = (. and a; = 0 , because Ji C (—cx),£] 

and I 2 C [ 0 , 00 ). If vi is continuous at a: = £, we have vi{tj = 0 , and (HTl) reduces to 

S( 6 ) = l + ^v[{^)v2{Q)lCot + 0{t). (43) 

Therefore, the leading order correction terms is at least of order The same holds for cases where 
V 2 is continuous at a: = 0. In this case, 

S( 6 ) = 1 - a;'(0) + Oi^). (44) 

For situation where Vi and V 2 are respectively continuous aX x = i. and a: = 0, the leading order 

correction term is at least of the order of In order to obtain the explicit form of this term, we 
must calculate Sq”'^ for n > 5. This is a tedious task. In the appendix we outline a different scheme 
to conduct this calculation. It gives the following improvement of (I4T1) . 

S( 6 ) = 1-^Ui(£)u2(0)/Co£^ 

~4^ {4^^iW^2(0)cri i[ui(£)u2(0) - u((£)u2(0)]/Co| 

{/c[ 6 ui(£)u 2 ( 0 ) - 4u((£)u2(0)]cri Vi{i)v2{0)a2 

+ i[- 4/c%i(£)u2(0) - v[{i)v2{0) + Vi{£)v2{0) + u"(£)u2(0) 

+ 4{ui(f')u2(0)^ -f- 'yi(^)^'y 2 ( 0 )}]/Co| £^ 

~ 32fc^ui(£)u2(0) - 6kv[{i)v2{0) 4fc{2ui(£)u2( 0 ) + v'l{£)v2{Q)} 

+ 16fcui(£)u2(0){2ui(£) U2(0 )}](Ti 4/c^ [5Ui(£)u2(0) - 2v[{i)v2{0)]CF2 

+ i[ 4 :k^{v[{i)v 2 { 0 ) - Wi(^)w2(0)}+ <(^)W2(0) -v[W2{0) 

+ Vi{i)v^2\^) - ^fH^)^^2(0) A{vi{£fv'2{^) - v[{£)v2{^f} 

- 14ui(£)u 2(0)K(£) - u'(0)}]/Co} t + 0{e). (45) 

In particular, whenever Ui(f') = ^ 2 ( 0 ) = 0, we have 
S(e)=l + ^u((£)u'( 0 )/Co^^ 

+ ^ {6fcu((£)u'(0) (Ji +i[v[{£)vm - <(£)u;(0)]/Co} t + 0{£^). (46) 


3 Application to Unidirectional Invisibility 


Consider the potential 


f 3 for 0 < X < L, 

\ 0 otherwise. 
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(47) 






where 3 , K and L are real parameters. This is the hrst-known example of a complex scattering po¬ 
tential that displays unidirectional reflectionlessness and invisibility [SI [HI UHl lUl 112 ] (for sufhciently 
small 3 .) In particular, if 

O'Trm 

(48) 


, K 2717X1 


for some positive integer m, then the reflection and transmission amplitudes of (147|) take the form 

P!: 


R‘ = 0(f), = +0(j’), r= l + m5"<=V + C>(3=), 


(49) 


where 


: = 


-iL"^ 




^( 2 ) .= 


iL^ 


Svri 


^( 2 ) _ 




dvrm^ ’ 327r3m^ ’ 4 

According to (09]), the potential (071) is left-reflectionless if we can neglect 0 ( 3 ^), and left-invisible 
if we can neglect 0(3^). Following (TU] we call these properties perturbative unidirectional reflec- 
tionlessness and invisibility, respectively. 

Let VK,m{x) denote the potential f07|) subject to the condition, 

4:77m 

L = Lm -.= —^, 

and consider the situation where 

vi{x) := Vk,i{x + Li- i), V2{x) := Vk,i{x), v{x) := Vi{x) + V2{x). (50) 

Clearly, the support Ij of Vj are given by 

h = [-L,+i,e], h = [0,Li], 

so that Ji n /2 = [0,f]. It is also easy to see that for i = 0, v{x) = Vk, 2 {x + Li). Because are 

invariant under space translations na, v{x) is also perturbatively reflectionless from the right for 
£ = 0. In the following, we use the results of the preceding section to explore the extent to which 
this property is violated for £ 7 ^ 0 . 

First, we compute the transfer matrices for Vj up to and including terms of order 3 ^ for 
k = K/2. We can easily do this by realizing that under a translation, x ^ x — d, the transfer 
matrix M of any given potential transforms according to M —)> g*fc(i(73_ Equivalently, we 

have R’^ RJ' and T ^ T, [12]. In view of (| 6 |), fl48|) . (09|) and fl50|) . these imply 


Ml = 
M 2 = 


1 + =^(2)^2 £-'^^(.^(1)3 + ^(2)^2) 
0 1 - ^( 2)^2 
1 + ^( 2)^2 ^( 1 )^^^( 2)^2 
0 1 - ^(2)^2 




C?(3^). 


(51) 

(52) 


Next, we compute S(e). It is easy to see that because vi{£) = ^ 2 ( 0 ) = 3 and v[{£) = ^ 2 ( 0 ) = *-^3; 
Eq. (0T]) gives 


S(£) = 1 -j|.K( 0 )(” 
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ai£^ + 0{£^). 


(53) 












Here we have also made use of 

Substituting fl^T]) - in ffTS]) . we find that for /c = i^/2 the transfer matrix of the potential v 
is given by 


M 


M 2 M 1 + S(e) - 1 + 0{f) 

' 1 + 2 (^( 2 ) ^ ^( 2))^2 ^ ^ 

1 _ 2 (^( 2 ) + ^( 2)^^2 


(54) 


where 


^( 2 ) _ 





±2ie + 
6 ^ 


(55) 


With the help of (E]) and fl541) . we can compute the reflection and transmission amplitudes of this 
potential for k = J’i'/2. These have the form 

= ^iV + (^(3^^') = |^He3 + e^)+0(3^6S), (56) 

RT = (1 + + C>(3^ t) 



T = 1 + V ^27r-^^ + 0 ( 33 , 6 ^). (58) 


As seen from these formulas, small overlaps between the support of vi and V 2 do not affect the 
unidirectional invisibility of their sum up to order 3 ^, but it slightly violates their unidirectional 
reflectionless up to order 3 ^. This is, in a sense, a sign of robustness of the phenomenon of unidi¬ 
rectional invisibility. 

In typical optical applications, v{x) = k‘^[l — n(a:)2], where n (x) is the refractive index of the 
medium [ 13 ]. This suggests that | 3 |/A '2 is of the same order of magnitude as | n(x)2 — 1|. Because 
for non-exotic material the real part of n is much larger than its imaginary part, the condition 
| 3 |/A '2 X corresponds to optical media with Re( n) ~ 1. To realize optical potentials of the form 
(j 37 j), one must dope a host medium with Re( n) ~ 1 and try to generate the needed loss/gain prohle 
by properly pumping it. This would necessarily involve various errors including those related to 
the location and size of the pumped region. Our results show that these errors do not obstruct the 
unidirectional invisibility and reflectionlessness of the sample. 


4 Concluding Remarks 

The composition or group property of transfer matrices is the main reason for their popularity 
and usefulness. This property applies whenever one wishes to determine the scattering features of 
a potential that is the sum of two or more constituent potentials with mutually disjoint support. 
In the present article we have explored a way of extending this property to a pair of constituent 
potentials Vj with overlapping support. 

For situations where the support of Vj intersect in a hnite interval [0, £], the standard composition 
rule, i.e., M = M 2 M 1 , generalizes to M = M 2 SM 1 , where S depends on i as well as the behavior 
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of Vj in [0,£]. Assuming that Vi are analytic functions in (0,-^), we can compute S in a power series 
in e = ki, where k is the wavenumber. We have shown that if Vi and V 2 are respectively analytic 
functions aX x = ^ and a: = 0 , so that Vi{t) = ^ 2 ( 0 ) = 0 , then the leading order term in the expansion 
of S — 1 is of order e®. Otherwise this term is of order or e^. 

Our results reveals another interesting fact regarding the scaling properties of the power series 
expansion of S. The coefficient matrices for n > 0, which determine this expansion, scale at 
least quadratically under the scaling of the potential: Vj —)■ avj, where a is a constant. This shows 
that the hrst-order perturbation theory is not capable of detecting the correction factor S — 1 , i.e., 
the standard composition rule for transfer matrices, namely ([7]), applies whenever the hrst-order 
perturbation theory is proven reliable. 
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Appendix: Power series for S(e) 

In this appendix we outline a method of evaluating the coefficient matrices appearing in the 
e-series expansion fl22|) of S(e). Because S(r) := M 2 (r, 0)“^M(r, 0)Mi(r, 0)“\ we should examine 
the e-series expansion for M(e, 0). It turns out that it is easier to work with U(r) := M(u 0), 

which satisfy 

S(e)=U 2 -'(e)U(e)Uri(e)e*^^T (59) 

Therefore, we hrst study the e-series expansion of U(e) and U(e)“^. 

Because 

/C(r) = e-*"^^/Coe'"^^ /Cq = (J3 + * (T2 = 

U(r) fulhls the matrix Schrodinger equation, 

iU(r) = H(r)U(r), (61) 

for the Hamiltonian 

H(r) := M;(r)/Co - 0 - 3 . (62) 

We also have U(0) = 1 . Therefore, U(e) = Texp [—i (irH(r)}. 

The Hamiltonian H(r) has been introduced in [10] and studied extensively in [TT] (See also [TT].j 
It enjoys the following useful property. 



U{rf = asiiiT) as. (63) 

This relation is a consequence of the fact that /Cg = as JCq as- As noted in [TU], /C(r) is a < 73 - 
pseudo-Hermitian matrix [T5], i.e., /C(r)l = ( 73 /C(r) 
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In view of fIM]) and fl62l) . we can easily verify that 


U(e) = (^lexpj Si-iH(r) e| ^ ^ 


n=0 


where 


cr := 


— i 


T = 0 




dr -i[w{T))C{oy - as] 


Similarly, we can show that 


f 

U(e)-i = V-C„ = (J2U(ef (J2, 

n\ 


Tl=0 


where 


C„ : = 


^ + iHfr') 


r=0 


= 0 - 20 ^ 0 - 2 . 


Here the last equality follows from fl63ll and fl65ll . 

In order to determine C„, we introduce the projection matrices: 


A : = 


1 1 
0 0 


r ; = 


1 -1 

0 0 


and note that, according to fl6Tl) . U(r) satishes 
*AU(r) = -rU(r), 


r=0 


'irU(T) = /(r)AU(T), 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


where /(r) := 2w{t) — 1. It is not difficult to see that rU(r) and AU(r) determine U(r) uniquely. 
Indeed, every 2x2 matrix A satishes 

i [(A A + r A) + a I (A A - TA)] = A. (70) 

We, therefore, proceed to explore AU(r) and rU(r). 

Clearly, AU(r) has the form 


AU(r) = 


u+{t) M_(r) 
0 0 


(71) 


for a pair of complex-valued functions m±, and in view of Eq. fl6^ and U(0) = 1, it satishes 

a2AU(r) = /(r)AU(r), (72) 


= ir. 


T=0 


(73) 


AU(0) = A, S,AU(t) 

These in turn imply that u±{t) are the solutions of the second order linear homogenous equation: 

u{r) - f{T)u{T) = (74) 

that fulhl the initial conditions: 


■*^±(0) = I5 ^±(0) = 


(75) 
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Next, we obtain a power series solution of fl7^ . To this end we express derivatives of u in the 
form 

MW(r) = p„(r) u{t) + g„(r) h(r), (76) 

where and Qn are a pair of auxiliary functions. In light of fl7i|) and fl76|) . they satisfy 


Pn+ii-r) 


where n > 1 and 'D[t) := 


= V{t) 


Pn{T) 

(Inir) 


’ Pli'r)' 


■ 0 ■ 

. <Pir) . 


_ 1 _ 


dr fir) 

1 dn- 


. In particular, 


relation in fl76]) and employing flTS]) . we find 


Pn+l{r) 

. Qn+iir) 


= V'^(t) 


0 

1 


M±(e) = 1 + 5^ 


dt 

(^ + 1 )!’ 


(77) 

. Using this 

(78) 


n=0 


where for all n > 0, 


h 1 ±i]v(T)’ 


(79) 


r=0 


Now, we substitute fl78l) in fITT]) to obtain AU. This together with fl6^ and flTOj) yield 


U(£) = 1 + 5 ^ 


^n +1 


* ^n+1 *^71 * ^n+1 


(80) 


^ 2(n + 1)! [ dn + i 

dn+l 

Comparing this relation to fIM)) . we hnd that the coefficient matrices C„ of Eq. flH^ have the form 



dn * '^n+1 *^71 * *^77+1 

< + ^4+1 + 


(81) 


This relation together with fl66|) and (l67)) allow us to compute the e-series expansion of U7^(e) and 
U7^(e). Using these and dSU]) in yields the desired expression for S(e). In the following we 
apply this method to compute for n < 6. 

First, we employ fl79|) to compute 


do =±i, df = f, df = f ± if, 

df = f + f±2tf, 4 = /(3) + 4//± ^(3/ + f), 

df = + 4/2 + 7 / / + ± z (4 + 6 / /), 

4 = /(') + 15 // + 11 //(=^)+ 9 / 2 /±^( 5 /(^) + 10/2 + 13 // + / 3 ). 


Here and in what follows / and its derivatives are evaluated at r 


0. Using these equations 
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in flHTl) . we obtain 


Ci = i 


a3--(/ + l)/C(0) 


C2 = /l-^/C(0), 


Cs = ^ /(31 - (Ti) - ^ [(/ +n m - mof 

C4 = (2/ + f)l-fa^ + tf }C{0f - t (2// + 

C 5 = ^ + 2 //) 1 - /(3) + / /) (Ti 

+ ^(/^'^ + /' + 4/^ + 7// + 3/ + / 2 ) (72 
+ f+ 4/^ + 7//-3/-/2) (73, 

C 6 = (3/(^) + 7/2 + 10 // + f)l- (2/(^)+3/2 + 3//) (7i 
+ ^ + 15 / / + 11 / + 9/ f + 4 /(3) + 6 / /) (72 

- ^ + 15 / / + 11 / +9ff-4 /(3) - 6 / /) ( 73 . 

Substituting these in fIMl) and fl671) . and using fl 66 l) . we can determine U(e), Ui(e)“^, and U 2 (e) 

and compare the result with fl 22 l) . we can compute Sq 


If we insert the resulting expressions in 
for n < 6. The result is 

sW = = 0, = -4zn;i(0)n;2(0)/C(0), 

Sg'^^ = — 16 t(;i(0)t(;2(0) (7i — 4i t(;i(0)t(;2(0) + 3 t(;2(0)t(;i(0) /C(0), 
16 t(;i(0)t(;2(0) + 6 t(;i(0)t(;2(0) (7 1 — 32 t(;i(0)t(;2(0) (J 2 


(n) 


(82) 

(83) 


S® = - 4 


q(6) _ 

Oq — 


- 4i 4 m;i(0)m;2(0) ( 2 m;i( 0 ) + 2 m; 2 ( 0 ) - 1 ) + m;i( 0 )m) 2 ( 0 ) 
+ 6m)i(0)m;2(0) + 4m;i(0)m;2(0) /C(0), 
16m;i(0)m;2(0){2m;i( 0) + ta2(0) - 1} + 15 m;i(0)m;2(0) 


(84) 


4t(;i(0)M)2(0) + 20 M)i( 0 )t(; 2 ( 0 ) 


(7i - 80 


Wi(0)t72(0) + 2 m;i(0)m;2(0) 


(J 2 


— Ai 


4t(;i(0)t(;2(0){17M;i(0) + 7^ 2 ( 0 )} + 4t(;i(0)M;2(0){2 m;i(0) - 1} 

+ 20m;i(0)m;2(0){2m;2(0) - 1} + 5m;i(0)m)2(0) + 10 {m)i(0)m; 2(0) + taf ^(0 )m;2(0)} 

+ '«;i(0)m;^^^(0) /C(0). (85) 

Next, we recall that we can assume, without loss of generality, that ai(e) coincides with its left 
limit at X = e. This allows for expanding tai(r) and its derivatives in Taylor series about r = e for 
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T < e. In particular, we have 

tyi(O) = wi{e) - wi{e) e + ^ wi{e) + 0(6“^), 

2 0 

m;i( 0) = M;i(e) - M)i(e) e + ^ + O(e^), 

m)i( 0) = M)i(e) - wf^(e) e + O(e^), M;f^(0) = wf^(e) + 0{e). 

Substituting these equations in (182|) - (l85l) . inserting the result in (l22l) . and noting that t = kx and 
Wj{T) = Vj{x)/2k‘^, we obtain fl45D . 
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